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ABS TRACT 


Several studies of rawinsonde component errors in actual flight 
tests are discussed. RMS relative errors between radiosondes are 
assumed for meaSured parameters. These errors are then used to derive 
RMS relative errors for other parameters such as wind speed and direction, 
height, virtual temperature, and mixing ratio. Expected RMS errors in 
measured and derived parameters are discussed in relation to the data 
obtained from the Oklahoma National Severe Storms Laboratory upper-air 
mesonetwork. It is found that many analyses and computations which are 
useful in synoptic-scale meteorology are not practical on the mesoscale. 


Wind data, when, processed by computer, yield the most valid results. 


DUDLEY KNOX LIBRARY 
Library : NAVAL POSTGRADUATE SCHOOL 
U.S. Naval Postgraduate School MONTEREY, CA 93943-5101 
Monterey, California 93940 


TABLE OF CONTENTS 


Te IM DROID TL A Go oe i nee” i 
de. EVALUATION OF RAWINSONDE COMPONENT ERRORS ......... 14 
A. PRESSURE ee ee, i es 14 

Be teenie... MM 5 ssi se ee a) we eel 116 

GC. REEVE eMUMEDITY a. © Se. . . eee. ww ww mm 18 

D. WIND MEASUREMENT SEM 05 5 oe race e020 are bie.” Rare 20 

ITT. RESU LTANESERR@ERSeIN DERIVED DATA . . . . «6 « » «© 1 «6 ew ww 26 
A. TEMPERATURE ERRORS a 26 
DMO MAIE mE eR IG RO Gers e 6 2 « 6 ws 6 6 6 eo 8 le 27 
Tweens Ratio andespeeture Humrdity . . 6. « « 6 6 27 

Pee Veet rem peraeUre (ssa. sate el wlth 6 8 we 28 

Se werelnt Temperature Wi . Wm. . ww wk lw lt lt lt wl 29 

C. HEIGHT OF STANDARD ISOBARIC SURFACES ee eee ee eS 30 

1. Error in Thickness Due to Temperature Error ... . 31 

pas Penola Thickness Due to Pressure Error « Den Eee 32 

3. Total Error in Height of Isobaric Surfaces ee ee 33 

D. HEIGHT OF THE RADIOSONDE See se is = 10) 2o8 <o. Rowe a Es 55 
MONDO | jes «ees se fs Se Oe we es we le we US 37 


IV. USEFULNESS OF VARIOUS ANALYSES AND COMPUTATIONS ON A 


OO UNAM Un. cyt snes © % 6 4 eee is 8 @ heel ew ww we EE 
POG EN IMONS hs Goes a ws 8 oe Wo © od eles oe Me eee we 45 
B. TEMPERATURE ANALYSIS ee he eet et ee has oS 
ree OM OMM tO 4. «Secs & ww 6 we ws Re se lw a9 
Dee Me oOsG ath WEND ACCURAGY . . 2. 2s 6 we ew se te 64 


E. ENERGY COMPUTATIONS 
F, ISENTROPIC TRAJECTORIES 
G. VORTICITY AND DIVERGENCE 


ve CONCLUSIONS | 


yi. SUGGESTIONS FOR IMPROVEMENT IN RAWINSONDE ACCURACY . 


BIBLUIOGRAPHY=:.. ... oie...) i cnoue 
INITIAL DISTRIBUTION LIST 


FORM DD 1473 


65 


66 


68 


ua 


74 


76 


v2 


81 


LIST OF TABLES 


Table Page 
I RMS Error in Dew Point for Sample Sounding 29 
II Standard Deviation in Layer Thickness Due a 


to Mean Virtual Temperature Error of 0.5C 


ied Standard Deviation in Height of Isobaric 32 
Surfaces Due to Mean Virtual Temperature 
Error of 0.5C 


IV Standard Deviation in Layer Thickness Due 33 
to a Pressure Error of 2.0 mb 


V Standard Deviation in Height of Isobaric 33 
Supiaces Due to a Pressure Error of 2.0 mb 


VI Standard Deviation in Height of Isobaric 34 
Surfaces (meters) including Effects of both 
Temperature and Pressure Errors 


VII Standard Deviation (meters) for Height of 36 
the Radiosonde at Selected levels 


VIII Standard Deviation of Position Error 40 
(meters) Due to GMD-1 Measurement Errors 


IX Standard Deviation of Position Error and 41 
Wind Speed due to GMD-1 Measurement Errors 
and Radiosonde Height Errors (meters) 





Figure 


2a 
2b 
Ze 
2d 
2e 
Da 
3b 
oc 
3d 
4a 
4b 


4c 


LIST OF ILLUSTRATIONS 


Position Error in Range Caused by Error in 
Elevation Angle of GMD-1 


850 mb Height Analysis 

500 mb Height Analysis 

300 mb Height Analysis 

Adjusted 850 mb Height Analysis 
Adjusted 300 mb Height Analysis 

850 mb Temperature Analysis 

500 mb Temperature Analysis 

300 mb Temperature Analysis 

Adjusted 500 mb Temperature Analysis 
850 mb Dew Point and Mixing Ratio Analysis 
500 mb Dew Point Analysis 


300 mb Dew Point Analysis 


Page 
38 


47 
48 
49 
51 
52 
54 
55 
56 
57 
60 
61 


62 





, 


: Pe 








}: «ee4eee8 T°” 





ATC 


m sec 


mb 


mc 


SE 


LIST OF SYMBOLS AND ABBREVIATIONS 


Advective temperature change 
Specific heat at constant pressure 
Horizontal displacement of radiosonde 
Vapor pressure 

Saturated vapor pressure 
Elevation angle 
Acceleration due to gravity 
Thickness (of a layer) 
Joules 

Latent heat of condensation 
Mixing ratio 

Meters per second 

Millibars 

Megacycle 

Saturated mixing ratio 
Montgomery stream function 
Nautical mile 

Pressure 

Specific humidity 

Balloon radius 

Dry gas constant 

Reynolds number 

Root mean square (see o) 


Static energy 


SR 


Slant range 

Time 

Temperature 

Dew-point temperature 
Virtual temperature 
Total energy 

Wind speed 

Geostrophic wind speed 
Height 

Potential virtual temperature 
Viscosity of air 
Density 


Standard deviation or RMS 


10 


I. ENTRODUCTION 


For several years the National Severe Storms Laboratory (NSSL) of 
ESSA has taken mesoscale observations and conducted research in an 
attempt to better understand the processes involved in severe local 
storms and therefore enable them to predict storm occurrence more 
accurately. In 1967 upper-air data were obtained by serial rawin- 
sonde ascents from each of nine stations in the southern Oklahoma 
mesonetwork on those days when severe storms were likely to occur in 
the area. GMD-I Rawin Sets were utilized to track the standard Weather 
Bureau radiosondes utilizing lithium chloride humidity elements. All 
soundings were terminated at 100 mb to insure a 90 minute interval 
between ascents. Data reduction was accomplished by computer utilizing 
essentially the same program as used for Project Stormy Spring. docu- 
mented by Kreitzberg and Brockman (1966). 

Using data provided by NSSL for three days in April and May, 1967, 
several graduate students at the Naval Postgraduate School have under- 
taken research projects to gain additional knowledge of atmospheric 
motions of this scale (cf. Anawalt, 1969; Van Sickle, 1969; Coleman, 
1969). All have encountered difficulty in analyzing the rawinsonde 
data. Analysis of mesoscale upper-air data is. newrelative to synoptic- 
scale analysis in meteorology, and there seems to be no proven tech- 
niques which apply to this scale. Attempts to apply synoptic-scale 
procedures have often failed or led to doubtful results. Ome cannot 
assume that the wind field is geostrophic or that balance equations, 


useful on a synoptic scale, will be valid. Also, gradients of moisture 


Ie: 


and temperature are expected to be highly nonlinear and much larger 
than is commonly measured using Synoptic-scale data. The main problem 
in analysis is: how much should the analyst smooth the data and still 
be reasonably certain of retaining significant features? 

At least part of the difficulty in analysis may be attributed to 
errors in rawinsonde measurements. In this case one is not so inter- 
ested in total error in temperature, pressure, and relative humidity as 
in the relative error or dispersion between two identical instruments as 
they would measure the same parameters. Since the upper-air mesonetwork 
as it existed in 1967 had as average distance between stations of 85 km 
(46 nm), the differences in both measured and derived parameters between 
stations sometimes result in unexpectedly large gradients. An attempt 
was made to discover what part of these gradients was attributable to 
real meteorological differences and what portion wee caused by relative 
errors between radiosondes. 

Sources of error and the results of several tests of radiosonde 
accuracy were studied and are discussed in Section II. Some of the 
factors which principally affect the total error are briefly mentioned. 
Scandaed deviations of temperature, pressure, and Petative huni 
between two similar radiosondes measuring these parameters at the same 


time and place are assumed. The standard deviation of relative error 


He 


is defined as: o= WT? where o is the standard deviation and alter- 
; ES , 
natively referred to as the root mean square (RMS) error, 2x is the 
sum of the squares of the differences, and N is the number of measure- 
ments of x. Resultant errors in derived data such as heights and 


mixing ratios are developed in Section III with the errors in wind 


data applicable to the synoptic network in the United States. Finally, 


r2 


these expected RMS errors in measured and derived parameters are dis- 
cussed in relation to the data obtained at the Oklahoma NSSL network. 
Computations and analyses which are useful on a synoptic scale are 
investigated for use on this mesoscale. 

The content of this thesis represents two main research objectives. 

They are: 

1. To develop an appreciation for the magnitudes of RMS relative 
errors in rawinsonde data as obtained from both the synoptic 
network and the NSSL mesonetwork. 

2. To discover what analyses and computations commonly applied to 


synoptic-scale data would be useful on this mesoscale. 


ie 


IIT. EVALUATION OF RAWLNSONDE COMPONENT ERRORS 


A. PRESSURE 

In the technical report, Accuracies of Radiosonde Data (1955), the 
Air Weather Service listed several sources of error and the tolerances 
for each in terms of standard deviation which were derived from flight 
Similitude tests. Sources of error itemized were standard barometer, 
Operating standard barometer, accuracy of pressure element, temperature 
effect, station barometer, and error in estimating fractional contacts. 
The standard deviations assumed for each of the above sources were 
statistically combined to give an expected total standard deviation of 
+ 2.12 mb up to 15 km (about 100 mb) for each radiosonde. A Signal 
Corps study of actual flight tests reported in the same publication 
resulted in a standard deviation of + 3.75 mb up to 15 km with smaller 
errors above that altitude. A value commonly ascribed to the standard 
deviation of pressure errors in a tropospheric sounding has been 
a = + 3.0mb. This value was used in this report by the Air Weather 
Service to compute height errors up to 100 mb. 

More recent investigations have shed some light on several possible 
sources of error. Hodge and Harmantas (1965) reported up to 3 mb 
difference for 16 standard Weather Bureau radiosondes between the origi- 
nal factory calibration and a recalibration at room temperature in the 
Weather Bureau laboratory. The recalibration of a like number of 
military AN/AMT-4B radiosondes revealed up to 6 mb difference. Most of 
the radiosondes tested, 95% of the Weather Bureau and 85% of the mili- 


tary, fell within + 2 mb. Some of these radiosondes were recalibrated 
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at a temperature 80 degrees colder than room temperature (about -55C) 
and the difference in pressure between warm and cold recalibration 
varied up to a maximum of 5 mb. Supposedly the baroswitch is tempera- 
ture compensated, but this is clearly inadequate. 

Since both the Weather Bureau and AN/AMT-4B radiosondes carry their 
batteries internally, this introduces a heat source to the other radio- 
sonde components. Andersen (1966) tested the heat buildup by placing 
a thermocouple one centimeter from the center of the interior of the 
baroswitch diaphram and measuring baroswitch and environment tempera- 
ture in a home freezer while a small electric fan prevented stratifi- 
cation. When these temperatures stabilized, the environment temperature 
was -28C and the baroswitch of the AN/AMT-4B was OC. This would imply 
that the baroswitch and, to some extent, the humidity element and 
thermistor would remain at some temperature which is warmer than the 
environment in a normal sounding of the troposphere. 

Corbeille (1966) compared the AN/AMT-4B with the AN/AMT-11 (a 403 
mc radiosonde used aboard Navy ships) in 50 flights of two radiosondes 
each. Spurious superadiabatic lapse rates found in his investigation 
were attributed to lag in the baroswitch of the AN/AMT-11. These lapse 
rates occurred between 400 mb and the tropopause over a time interval 
of one or two minutes. The lapse rates were not verified by the 
AN/AMT-4B trace even though the pressure did not decline at a steady 
rate. However, the AN/AMT-11 uses an externally mounted battery while 
the AN/AMT-4B, like the Weather Bureau radiosonde, uses an internally 
mounted battery. The theory advanced is that the heat from the battery 
compartment of the AN/AMT-4B which keeps the baroswitch much warmer 


than that of the AN/AMT-11 may cause errors due to temperature but at 
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the same time lessens the lag due to friction of the contact arm on the 
commutator bar and prolongs the relatively smooth functioning of the 
baroswitch. 

The pressure difference between the AN/AMT-4B and the AN/AMT-11 was 
computed each minute of 15 flights by Corbeille (1966). The AN/AMT-4B 
had a mean pressure 2.0 mb lower than the AN/AMT-11C, and a standard 
deviation from the mean of 0.9 mb resulted. The mean differences were 
evaluated at three minute intervals and showed a range of -3.5 to + 2.0 
mb. Hodge and Harmantas (1965) computed the root mean square of the 
pressure differences taken at one minute intervals for 16 flights of the 
Weather Bureau and AN/AMT-4B radiosondes. This standard deviation was 
2.l@mnb fommall 16 Elights andgsanged from 1.1 to 3.1 mb formeach flight 
The maximum pressure difference of 7.0 mb occurred twice in the lower 


troposphere. 


B. TEMPERATURE 

In the publication Accuracies of Radiosonde Data (1955) the Air 
Weather Service synthesized four studies of temperature errors including 
radiosonde comparisons by New York University and the Oklahoma City 
tests by the Weather Bureau in 1951. Assuming that errors were normally 
distributed around zero, the standard deviation of temperature error was 
computed at standard pressure levels and ranged from 0.65C to 1.3C in 
the troposphere for the NYU study. The Oklahoma City tests showed a 
standard deviation of 0.82C for all standard pressure surfaces from the 
surface up to 400 mb and 0.99C from 300 mb to 100 mb. It should be 
noted that this data represents the errors to be expected in a one 


instrument sounding, and therefore two instruments of the same kind 


16 


flown together would have a standard deviation of their differences of 


1.414 times the standard deviation of a one-instrument sounding according 


a 0 av oe Also, the data includes errors in pressure pro- 


Coro 
portional to the baric lapse rate. This error source will be investigated 
more fully in Section III. 

It has been common practice by analysts to assume a standard devi- 
ation of 1C for upper-air temperature data on isobaric charts. This is 
close to the results of all but the most recent studies. Corbeille 
(1966) in comparing the AN/AMT-4B and the AN/AMT-1IC found that in some 
soundings temperature differences greater than 5C existed, but the 
majority of soundings contained differences of less than 3C. The 
AN/AMT-4B was biased toward warmer temperature with the mean value of 
the differences being 0.6C and a standard deviation from that mean of 
0.9C. These figures should not be strictly applied to the Weather 
Bureau radiosonde because, as mentioned before, the AN/AMT-4B, like 
the Weather Bureau radiosonde, has an internal heat source while the 
AN/AMT-11C does not. Also slight difference in thermistors existed 
between the two radiosondes. However, the results are compatible with 
an assumed standard deviation of 1C for each instrument. 

Air Weather Service Pamphlet 105-3 (1967) specifies a RMS accuracy 
of 0.7C for the temperature element in the AN/AMT-4, the Weather Bureau 
radiosonde and the AN/AMT-12 (which contains a hypsometer and is some- 
times used by the U. S. Air Force). This would indicate an expected 
RMS error of about 1C between two radiosondes on the same flight; but 
this value applies up to 36 km which is considerably higher than the 
tropopause, the limit of interest for this study. Daniels (1968), 


Reynolds and Lamberth (1966), Ney et. al (1961), and Badgley (1957) 


fey 


indicate that radiation errors, lag constants, ventilation rates, and 
configuration effects contribute more to temperature measurement 
errors with an increase in altitude and considerably affect the total 
accuracy. Since relative accuracy is of interest here, one would expect 
somewhat less error if soundings were limited to 100 mb (about 15 km). 
Hodge and Harmantas (1965) computed the RMS differences for each 
minute of 16 flights which resulted in a RMS value of 0.51C for all 
flights and a range of 0.38C to 0.8/C for each sounding. The tempera- 
ture differences at the mandatory pressure levels were also computed 
with 53% of the daytime differences and 52% of the nighttime differ- 
ences within + 0.5C. These figures result in a standard deviation of 
0.7C for the standard isobaric levels. The standard deviations of 
0.5C at the same time and 0./7C at the same isobaric level were con- 
siderably better than one would expect which indicates that differ- 
ences between compatible instruments are much less than the total 


error ,of thewimst rument.« 


C. RELATIVE HUMIDITY 

The lithium chloride humidity element was being phased out in 
favor of the carbon element, but it was the LiCl element which was 
used in the radiosondes at the Oklahoma NSSL network in 1967. Ma jor 
shortcomings ascribed to the LiCl element are insensitivity, large lag 
constants, and the reluctance to meaSure high humidities in the presence 
of rain or cloud droplets. O'Conner (1952) indicated that relative 
humidity is measured with respect to water; and, therefore, lower 
readings should be expected in cirrus clouds. In that case 70-80% 


relative humidity may be accurate. For example, at -20C 83% RH with 
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respect to water is equivalent to 100% RH with respect to ice. Also, 
condensation on the element leads to error at the time of occurrence 
and during the rest of the sounding. Mathews (1964) ascribes the low 
readings in clouds to an exothermic liberation of heat in the lithium 
chloride film as it absorbs moisture. 

Mathews (1964) summarized the development of the lithium chloride 
hygrometer and presented the testing specifications required of the 
manufacturer. By 1962 all elements were tested at room temperature to 
+ 3% of a reference relative humidity. Four out of 100 were tested 
against the standard calibration between 20% and 95% ascending and 
descending relative humidity at eight points. Errors were not to 
exceed 7% RH at one point and 5% at all other points. One out of 100 
elements were tested for polarization by switching it on for 15 seconds 
and then off for 30 seconds for a total of 30 minutes at a constant 
relative humidity. Maximum drift allowable was 6% relative humidity. 
One element out of 1,000 was tested for time lag by decreasing the 
relative humidity from 90% to 40% and allowing a maximum of 15 seconds 
for 90% of the change with a wind speed of 250 meters per minute (the 
approximate ascent rate of the standard radiosonde). The. first three 
test specifications indicate a standard deviation of those errors com- 
bined to be about 3% relative humidity at room temperature. 

Large lag constants would cause large errors in relative humidity 
especially if the relative humidity were changing rapidly. Mathews 
(1964) also reported a study of low temperature lag characteristics 
done by Jones and Wexler (1960). At a wind speed of 160 meters per 
minute the response time to accomplish a 63% (1 - 1/e) change in 


relative humidity ranged from 51 to 74 seconds at -20C and 120 to 480 
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seconds at -40C. Ference (1951) reported response times of four seconds 
at -25C, 15 seconds at OC, and 120 seconds at -30C. This implies that 
at 500 mb and higher this element may indicate relative humidities 
which were sensed hundreds of meters below. For computations at standard 
pressure levels or even averaging over small layers the lag of the lithi- 
um chloride hygrometer may cause relatively large errors. 

It has been standard practice by analysts to assume a standard 
deviation of + 5% RH for temperatures above OC and + 10% between OC 
and -40C. The corresponding standard deviation of the difference 
between two instruments on the Same sounding would therefore be 7% and 
14% respectively. Hodge and Harmantas (1965), using the carbon element 
in the AN/AMT-4B and the lithium chloride element in the Weather Bureau 
radiosonde, found that 79% of the differences at each minute were with- 
jaa + 10% RH. This figure converted to a normal curve results ina 
standard deviation of 8% RH. Corbeille (1966) reported essentially the 
Same overall results with the errors increasing with decreasing pressure. 
Unfortunately there have been no tests in the last decade with two 
lithium chloride hygrometers in the same flight. Therefore a standard 
deviation between two sensors of 5% RH above OC and 10% RH below OC 


Witlilebenacaumecdes 


D. WIND MEASUREMENT 

Wind is commonly measured by tracking the radiosonde in its ascent 
by a GMD-1 Rawin Set and recording azimuth and elevation angles at set 
time intervals. The height of the radiosonde is computed from the 
integration of the hypsometric equation; and by using azimuth and 


elevation angles, the position of the instrument over the curved earth 


yA 


is thus determined. By comparing two positions at either end of an 
elapsed time one arrives at a mean wind speed and direction for the 
layer represented by the heights of the radiosonde at either end of 

the time interval. GMD-1 wind data are not reported when the elevation 
angle is less than six degrees above the horizon or any prominent object 
or when the azimuth angle is within six degrees of any prominent object 
on the horizon because of reflected or refracted signals. The RMS 
accuracy of GMD-1 elevation and azimuth angles is thought to be + 0.5 
degrees (cf. Danielsen and Duquet, 1966). 

The accuracy of the wind obtained from such measurements or from 
measurements of any free balloon is directly related to the validity 
of two rather gross assumptions. The first assumption is that the 
balloon responds completely to the action of the wind. That is, no 
relative motion exists between the balloon and the air mass except that 
provided by the balloon's ascent. Consider what happens to the balloon 
as it encounters changing winds in its ascent. Naturally the balloon 
trains' mass, momentum, and aerodynamic drag and the air density and 
vertical shear of the wind are factors which must be considered in 
determining how fast the balloon train does attain a negligible hori- 
zontal relative motion. Near a jet stream the vertical shear is large 
enough to cause some wind errors, but elsewhere the errors of this type 
could be considered negligible. 

In the last decade there has been increased interest in small scale 
fluctuations of the wind because of its importance in launching rockets 
and space probes. As a result there have been many studies of balloon 
motion relative to the air and another type of error has been brought 


to light. There exists a self-induced balloon motion which has been 
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described as both a regular spiral or perturbation vector about a mean 
wind vector and highly erratic depending upon Reynolds number (cf. 
MacCready, 1965). Reynolds number is defined as Re = 2pVr/u, where op 

is the density of air, UW is the viscosity, V is the speed of air rela- 
tive to the balloon, and r is the balloon radius Causauae? L965)e4 With 
Supercritical Reynolds number the motion is erratic while subcritical 

Re permits more regular balloon oscillations. 

Murrow and Henry (1965) investigated the behavior of several types 
and sizes of balloons including the standard radiosonde balloon in still 
air inside a 220 foot high blimp hanger at Lakehurst, New Jersey. They 
concluded that horizontal velocity was about one-half the vertical 
velocity, and that large roughness elements added to a balloon surface 
resulted in reduced amplitudes of horizontal motion and reduced wave- 
lengths of this motion in the vertical. Scoggins (1965) investigated 
the behavior of ROSE and JIMSPHERE balloons and concluded that increasing 
the number of large roughness elements on the JIMSPHERE reduced the 
oscillations. McVehil, Pilié, and Zigrossi (1965) used doppler radar 
to confirm that the type of oscillations of one- and two-meter 
spherical balloons are dependent upon Reynolds number. MacCready 
(1965) summarized previous work. He suggests the use of the JIMSPHERE 
to define small scale variations in the wind field from the surface to 
18-20 km. The standard radiosonde should go through a critical Reynolds 
number regime between 8 and 12 km to a subcritical regime where oscil- 
lations become more regular. 

Lhermitte (1967) utilized doppler radar in his study of balloon 
motion. In 45 flights he used neoprene balloons of diameters from 0.3 


to about two meters. His conclusions imply that a radiosonde balloon 
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would develop erratic horizontal displacements of as much as six or eight 
meters in the supercritical Reynolds number regime (below 8 km according 
to MacCready (1965)). At higher altitudes in the subcritical regime 
balloon displacements were regular and about the same amplitude as the 
balloon diameter (about 3 meters) and the vertical wavelength would be 
about 50 meters. Therefore, the winds derived from radiosonde meaSure- 
ments averaged over two minutes (about 600 meters in height) are rela- 
tively unaffected by the horizontal motions, but they only reveal the 
gross features of the wind field. As the averaging interval decreases 
these self-induced horizontal motions take on added significance. 

The second assumption is that the position of the balloon or radio- 
sonde can be accurately determined. This results in the most error in 
routine rawinsonde wind reports. The Air Weather Service and Sandia 
Corporation conducted a series of wind tests in aaity Pose at Salton 
Sea test base. These tests incorporated wind data from the surface to 
at least 100 mb. Smith (1954) reported the results of those tests. 
Winds from three askania phototheodolites were used as a standard in 
six radiosonde flights to test the accuracy of two SCR-584 radars and 
two GMD-1(A)s. The Askanias' reportedly have a position error of less 
than three meters at 12 km. Therefore the use of Askania-derived winds 
as a standard appears logical. Although balloon positions and elevation 
angles were not compensated for curvature of the earth, the maximum 
error incurred was reported to be about one m — Balloons were 
underinflated in an attempt to test accuracies at the lowest elevation 
angles, but none were read at less than 10 degrees. The standard 
deviation for wind speed as measured by the two GMD-1(A)s was 2.2 


m Reet for all flights and all altitudes and 3.4 m ger > at altitudes 
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above 7.5 km. Standard deviations from all stations were directly 
correlated with height, slant range, and wind speed. 

Danielsen and Duquet (1966) developed and tested a machine method 
for computing winds from GMD-1 tracking and radiosonde data. FPS-16 
precision radar derived winds were used as a standard. The RMS 
accuracies for the radar are almost one order of magnitude smaller 
(0.01 degree in elevation and azimuth and five yards in slant range) 
than GMD-1l and radiosonde accuracies. Radar readouts were obtained 10 
times a second to determine small scale variations in elevation and 
azimuth. Winds were computed over 2, 10, 20, and 40 second averages. 
The power spectra of the winds revealed a large en peak at a 
period of nine seconds which was attributed to the natural period for 
self-induced balloon oscillation. The two second averaged winds 
showed fairly regular oscillations of + 2m as It is interesting 
to note that the authors had to use 40 second wind averages to eliminate 
these high frequency oscillations from the 2-second averages, and the 
40-second averages were used in later comparison with GMD-1 and radio- 
sonde derived winds. 

GMD-1 readings were taken every six seconds which is the maximum 
rate. Deviations from the mean elevation angle appeared to be 0.2 
degrees or four times the commonly accepted standard deviation. To 
Suppress this error the authors suggested that an average over at 
least six points (36 seconds) was necessary. il meaner the wind was 
computed over an interval equal to the barostat contact interval. 

This required averaging five or six azimuth and elevation angles near 
the surface and 10 or 12 at 100 mb. Increasing the number of obser- 


vations in each average helped to compensate for decreasing accuracy 
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in height and elevation angles. By using an average elevation angle 

and height, the distance over the curved earth was computed. With this 
distance and average azimuth angle the position was obtained. Then by 
taking finite centered differences of successive mid-point positions the 
winds were calculated at each pressure contact. The time interval varied 
from about 30 seconds to more than 60 seconds. 

In one ascent with fairly light winds the GMD-1/radiosonde derived 
winds were within one m seen and two degrees of the FPS-16 winds. For 
another sounding with a strong (68 m eee +) jet stream the GMD-1/radio- 
sonde winds accurately described the wind field at lower levels; but 
above the jet core (13 km) and at low elevation angles, errors of three 
or four m ie were evident. Two ascents resulted in elevation angles 
of less than 10 degrees. Maximum error in elevation angle was 0.7 
degrees or 1406 making it necessary to fit a fourth order polynomial 
to smooth the elevation angle data in order to compute reasonable 
winds. The authors concluded that GMD-1 derived winds closely re- 
sembled FPS-16 winds when GMD-1 angles were measured 10 times a minute, 
carefully read and processed by computer, and elevation angles remained 


more than 10 degrees above the effective horizon. 
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III. RESULTANT ERRORS IN DERIVED DATA 


It is necessary to assign some RMS error values to relative humidity, 


temperature, and pressure in order to compute errors in derived data. 


Again it must 


be emphasized that these values are not indicative of 


absolute accuracies but of relative error or the RMS dispersion between 


Weather Bureau radiosondes measuring the same parameters at the same 


time and place. 


Hodge and 


Harmantas (1965) concluded that the AN/AMT-4B and the 


Weather Bureau radiosondes were indeed compatible because differences 


noted between 
expected from 
facturers for 
utilizing two 
clusion takes 
consideration 
the following 


tionw@: 


Temperature: 


Pressure: 


Relative Humidity: 


these two radiosondes were no greater than could be 
procurements for different years or from different manu- 
the same agency. Since there have been no recent tests 
Weather Bureau radiosondes on each sounding, this con- 

on added significance. One must give their results great 
in assigning values for further computation. Accordingly, 


RMS relative errors will be utilized in future computa- 
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) 
oe 

7 
= 
o 


Ij 

V 

© 

a) 

|. 

oa 
“eA 


when clear of clouds 


=40C5< Te 9) GC, + 10% and precipitation 


A. TEMPERATURE ERRORS 


To discover temperature errors at some standard pressure levels, 


One must consider the pressure error and the prevailing lapse rate. 


The additional error at a particular pressure level will be: 
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dT -< Oo , where a = baric lapse rate at the level in question and 
4 cP E 
* = RMS error in pressure (2.0 mb). Using a sounding which is repre- 


sentative of conditions in the NSSL mesonetwork at 1830Z 28 May 1967 
results in additional errors due to pressure of from 0.17C at 850 mb 
to 0.46C at 200 mb. Since the lapse rate near 100 mb was essentially 
zero, there would be no additional correction for the 100 mb pressure 


level. The total temperature error on a constant pressure surface or 


at any specified pressure would be Op = /(o.5)° a0 at)” Hodge and 


Pp 


Harmantas (1965) computed RMS temperature errors at standard pressure 
levels and found the average to be 0./C. In the lower troposphere the 
error is probably less while at 200 mb the error is probably more than 


Om. 


B. MOISTURE PARAMETERS 


1. Mixing Ratio and Specific Humidity 


Since q = there exists a small relative error in assuming 


ek 
1+m’ 
that q =m. This error is about 1% for a q or m of 10 gm wen 
Since RH = ie an error of 5% in RH implies an error in m of 
% . 
0.05 m.. For example, consider a relative humidity of 50%, tempera- 
S 
ture of 23C and pressure of 1000 mb corresponding to a m. or V6.0 
gm een A 5% error in relative humidity corresponds to an error in 
mixing ratio of 0.9 em foie and therefore the actual mixing ratio 
probably lies between 8.1 and 9.9 em ert and the corres ponding dew 
pommteirom 10.5 to 13./C. At a relative humidity of 50%, a RMS error 
of 5% RH results in an error of 10% of the reported mixing ratio; and 


at 25% RH for subzero temperatures, a RMS error of 10% RH results in an 


error of 40% of the Reperted mixing ratio. Thus the error in mixing 
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ratio or specific humidity is tied to the saturation value - not the 
reported value. 
2. Virtual Temperature 
Virtual temperature is defined as; 
T. = i({)l + Cl org). 
where T is air temperature in degrees Kelvin. If mixing ratio is used 
in place of specific humidity, virtual temperature can then be written 


as 


T, = T+ 0.6l1mT. (1) 


One finds T = 285K and m= 7 gm on near 7/50 mb on the sample sounding. 
Therefore Ty becomes 286.22K or 1.22C higher than the reported tempera- 
ture. Taking the differential of equation (1) for error analysis: 

dT = dT + 0.61(mdT + Tdm). 
If one lets the differentials assume the values of the appropriate 
standard deviations, then dT = 0.5C =o 


T 
Me. z 12 gm erie Therefore mdT = 0.0035C and Tdm = 0.171C. Since 


and dm = 0.05m_ = 0.0006 at 


mdT < < Tdm, that term can be dropped. Even 0.61 Tdm is small (0.104C). 

If one assumes that dT = Om > then equation (1) now reads 

V 

On = 0.5 + 0.104 = 0.604C. However, Or is an average over the whole 
V 


sounding, and in the lower troposphere the dispersion in temperatures 
is less than average (Hodge and Harmantas, 1965). In the lowest layers 


of the tropical atmosphere 0.61Tdm may be large enough to increase On 
Vv 


Significantly, but in mid-latitudes the reduced dispersion in tempera- 
ture in the lower layers will compensate for errors in mixing ratio in 
determining virtual temperature. Consequently, in further computations 


O = 0.5C will be assumed, 


al 
V 
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3. Dew-Point Temperature 


T the dew point, depends upon the temperature, pressure, and 


d? 
relative humidity (or mixing ratio) and can be read from the thermo- 
dynamic diagram. Mitchell (1967) converted relative humidity errors 
of the same magnitude assumed here to dew-point errors for each two 
degrees of temperature from +40C to -40C in increments of 5% relative 
humidity. Also presented are two nomograms for the temperature ranges 
40C above and below freezing based upon interpolation of the tabulated 


data. In computing dew points, Mitchell assumed that RH a 2S, and 


S S 
since em the saturation vapor pressure, is dependent only upon 


a 
@ 


temperature, his results should be applicable at all pressures. 
Using the same sample sounding and Mitchell's tables, RMS 
errors in dew point are listed in Table I for selected levels. 
TABLE I 


RMS Error in Dew Point for Sample Sounding 





Pressure im (ep RMS Error in TCC) 
950 26.4 mero a we 
850 Pre. Pr a 
700 6.0 -2.9 aL 
500 -10.5 -24,2 ri a 
300 36 © => 99 + 4.3 


The errors are greatest for high temperatures and low relative 
humidity in each temperature range. The temperature error would cause 


an error in m. and therefore m and Ty at a particular pressure level 
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apart from the error caused by relative humidity. For example, if the 
temperature at /00 mb in Table I was in error by 0.7C and the relative 
humidity remained constant, the reported dew point would be about -3.5C 
or -2.5C corresponding to a temperature of 5.3C and 6.7C respectively. 
However, the RMS error in dew point due to errors in relative humidity 
would remain about the same (+ 1.3C), and the combined RMS error would 
be about + 1.4C. Therefore the temperature error is of little conse- 
quence in determining dew point accuracies and may be considered 
negligible. 

C. HEIGHT OF STANDARD ISOBARIC SURFACES 

One must Freer e tte between the heights of a constant-pressure 
surface and the height of the radiosonde at the time it indicated that 
pressure. While both heights are computed in the same manner, the 
errors in the latter heights are much greater and will be dealt with 
in the next section, 

The heights of constant-pressure surfaces are computed from the 
hypsometric equation using a mean virtual temperature for the layer. 
The thickness of each layer is added to the sum of the thicknesses of 
the layers below to arrive at a height for the top of the layer. The 
height of the first standard pressure level above mean sea level is 
computed by the barometric equation with less accuracy but of smaller 
thickness than the next layer, so that in computing height errors, it 
is assumed that the 1000 mb height is known exactly. The hypsometric 
equation which follows is derived from the equation of state for moist 


air and the hydrostatic equation (Haltiner and Martin, 1957): 


2 p 

R R 1 
-_ — ad — ee — a ~ 2 
Zour 2, : Jt, 6(1n p) Bi In - (2) 
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The subscripts 1 and 2 represent the bottom and top of the layer 
respectively and Et is the mean virtual temperature in the layer. By 
assuming that R and g are constants then the error in thickness of a 
layer bounded by Py and Po will be dependent upon the error in r., the 
baric lapse rate, z, assumed constant in the layer, and the standard 
deviation in pressure measurement, oT The error in thickness, dh, 


from Accuracies of Radiosonde Data (1955) is given as 


p 

R 1 Or 

dh == 1 (at : > 3 
Ears Sot OPP .” 


i. “Srrer ime Bhirckress Due to Temperature Prror 


The temperature component of the thickness error is 


Pp») — 
In— dT. (4) 
p> = 


dh,, = 


09 | 


Since one does not use the virtual temperature at any parti- 
cular pressure but only the mean qT in the layer, a standard deviation 
of 0.5C applies; and the results are shown in Table II. 

TABLE II 


Standard Deviation in Layer Thickness 
Due to Mean Virtual Temperature Error of 0.5C 


Layer (mb) dh,, (meters ) 


T 





ol 


If it is assumed that the temperature error maintains the same 
sign throughout the sounding, the standard deviation in the height of 
each pressure surface is obtained by adding the thickness errors. Table 
III gives the results of height error, Our? of several pressure surfaces 
with only the error in temperature considered. 

TABLE III 


Standard Deviation in Height of Isobaric Surfaces 
Due to Mean Virtual Temperature Error of 0.5C 


Pressure (mb) (meters ) 


Cnr 





2. Error in Thickness Due to Pressure Error 


The portion of the thickness error due to pressure is 


p 
dh = Es ein 1 of 


— Sie (5) 
p fc \peaepaae 
Because the sign of the error is not important in this application, a 
positive sign is assumed. Notice, however, that if the sign of the 
lapse rate changes with an otherwise constant pressure error, that 

. e. wOl 
the cumulative errors due to pressure would decrease. Using aa from 


the sample sounding and o, = + 2 mb, Table IV was compiled. In the 


sample sounding the tropopause occurred at about 150 mb with the 


SY 


temperature at 200 mb and 100 mb being equal. Since = is assumed 
constant in the layer, there is no error in the 200-100 mb layer; and 
——« additive effect of pressure error on height, in Table V shows 


the same value for 100 mb at it does for 200 mb. 


TABLE IV TABLE V 
Standard Deviation in Layer Standard Deviation in Height 
Thickness Due to a Pressure of Isobaric Surfaces 
Error’of 2.0 mb DWe to a Pressure Erre@t of 240 mb 


Layer (mb) ares) Pressure (mb) a 





oS. Sietal Error in Height of Bsobaric Surfaces 


Assuming that temperature and pressure effects are independent 
of each other and normally distributed Op» the total error in height of 


/ 2 2 
a constant-pressure surface in meters,’ is given by 7 (oO) aD ; 


and o from’ tables “LLL and V are included in Table VI for 


Ont hp 


comparison. 

Total error in height at standard constant-pressure surfaces 
is therefore mostly dependent upon mean virtual temperature errors. 
The assumption that the temperature and pressure effects on height 


errors are normally distributed and independent of one another seems 
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more plausible when one is considering relative errors between two 
radiosondes than when one is trying to obtain absolute errors. 
TABLE VI 
Standard Deviation in Height of Isobaric Surfaces 


(meters) Including Effects of 
Both Temperature and Pressure Errors 





Pressure (mb)| Our | OL 
1000 0 | 0 0 
oe | oe | 5.5 
500 a 610.1 : 3.) ll L0G 
300 17.6 ! | oe | 19.1 
200 23.8 Pe 12 Ome 26.8 
100 i 33.7 i 1279) Oeil 36.1 





The true height error in geopotential”meters is intluenced” by 
several factors which do not affect relative error. The fact that g 
is considered a constant in the hypsometric equation decreases the 
height of the 100 mb surface by about 40 meters. However, Ene Lact 
that the emissivity of the thermistor is near 0.3 on daytime flights 
probably compensates for the error in g by providing a bias toward 
warmer than true temperatures (Daniels, 1968; Corbeille, 1966). Also 
there is a small lag in thermistor response which could increase the 
temperature beyond the steady-state error (Badgley, 1957). At night 
the radiation error would cause the thermistor to record lower than 
ambient temperatures thus increasing the error in true height (Daniels, 


1968). While these are some factors which affect the true accuracy to 
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a considerable extent, one must assume that they are negligible in com- 
paring two radiosondes flown at the same time of day in the same geo- 
graphical area. One must assume that both thermistors are exposed to 
equal amounts of solar radiation. Otherwise the radiation error would 
influence the relative error in temperature between radiosondes. Of 
more interest are the random errors caused by deviation from standard 
in manufacture, lock-in errors during the baseline check, proper align- 
ment and readability of the recording equipment, and accuracy of data 
reduction. These types of errors seemingly are more apt to cause a 
normal distribution of errors in derived data, because any bias caused 


by inadequate design is common to all soundings. 


D. HEIGHT OF THE RADIOSONDE 

It was shown in Section C that errors in height of the constant- 
pressure surfaces were caused mostly by temperature errors and only 
slightly increased by the pressure error coupled with the baric lapse 
rate. However, if one tries to determine the height of the radio- 
sonde when it indicated a given pressure, then one must consider the 
error in meters represented by o. G-2emb) sat Eheesurfiace Ginequvestion. 
For instance, 2 mb in pressure at 500 mb represents 31 meters. However, 
for the applications considered in this section, one is not interested 
in relative error between radiosondes but rather that portion of rela- 

O 

time error which is attributable to a single radiosonde or —-, From 


J2 
the hydrostatic equation and the equation of state for moist air 


R dip 0 
dz = - —~“—. Replacing dz and dp by O__— and aes respectively, this 
& P 2p ge 
becomes 
i (+ /2) BB 
anil? IO eee Mat 29.29 “Ty (6) 
2p sa Pp a : = 
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Using the values for qT. from the sample sounding at the corres- 
ponding pressure, the height error corresponding to a 2 mb pressure 
error is obtained. Since wind and significant temperature levels are 
attributed to the heights of the radiosonde at the time they were mea- 
sured, this error in height due to error in pressure is frequently 
encountered. For a single sounding the RMS error in height due to 


pressure, O, , and the RMS error in height of a constant-pressure sur- 
Pp 


0 
face, —-, are combined to get standard deviation in radiosonde height, 





2 2 
(9,)° + (o,.) 
CDC Pp 6, «©The results follow in Table VII where 





RMS errors in some standard isobaric levels are used for OL: 


TABLE VII 


Standard Deviation (meters) for Height 
of the Radiosonde at Selected Levels 


Pressure (mb ) | 





As can be readily observed from Table VII, the error in pressure 
is the predominant error in height of the radiosonde. Although 
stratospheric measurements are beyond the scope of this thesis, it 
is interesting to note that above 100 mb the pressure error is probably 


smaller using the conventional barostat; but the error in height due to 
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presswre error, eae at 10 mb would be about five times as great as 

the error at 50 mb. Dowski (1961) showed that barostat-derived heights 
were fairly accurate in the troposphere with errors commonly about 20 
or 30 meters and a maximum error in six flights of 75 meters. However, 
at an altitude of 30 km, errors as great as 1 km were observed while 
hypsometer-derived heights at that altitude were generally less than 


200) meters in error. 


E. WIND ACCURACY 

Suppose one makes some rather gross assumptions in order to make 
some simple calculations and gain a feeling for the errors in radio- 
sonde position caused by uncertainty in azimuth and elevation angles 
of the GMD-l. First, assume that the height, z, is known exactly and 
that the standard deviation of the azimuth and elevation angles is 0.05 
degrees. Second, assume that the earth is flat; and, therefore, the 
effective elevation angle is also the true elevation angle. The rep- 
resentation of a spherical surface by a flat surface results in an 
error which increases with distance of the radiosonde from the tracking 
site. For instance, if the radiosonde were displaced 111 km (60 nm), 
the angle subtended by radii from the earths" center to the tracking 
site and radiosonde would be one degree causing the elevation angle to 
read one half of one degree too low. However, elevation angles to be 
used in this calculation will be 14 degrees and 26.6 degrees with a 
maximum displacement of 60 km. Therefore, the maximum error in eleva- 
tion angle is 0.27 degrees. 

In Figure 1 the following definitions apply: 


D = horizontal displacement 


a7 





Figure 1: Position Error in Range due to RMS Error in 
Elevation Angle of GMD-1. 
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SR = slant range 
E = elevation angle 
P = position of the radiosonde 


Z = height of the radiosonde above the GMD-1 


O, = standard deviation of elevation angle = 0.05° = 0.00087 
radians 
Ope = error in range on a horizontal plane due to RMS error in 


elevation angle 


Since Z, E, and O, are known, one first must find SR. This is 


given by Z/sin E. Therefore the distance PQ is given by Z0,/sin E. 


However, Opp? the error in range due to On is PQ/sin E. Therefore Ope 
is 
ZO 
E 
gS ——— (7) 
Z 
(sin E) 


Similarly the error normal to the line of sight is SRO, and is 


ZO 
Z A 
oa a (8) 


where On» 


-l1 
the ascent rate of the radiosonde is constant at five m sec, Further, 


the RMS error in azimuth angle, equals Op: Now assume that 


consider two cases where the mean wind speed displacing the radiosonde 


iswhOsandeZ20 m ociae so that tanE is 0.5 and 0.25 respectively. Opp? 


and o the standard deviation of position error due to GMD-l 


Oy? XG? 


measurements, are given for both cases in Table VIII where 


Oya = om + (04). 


Since wind is computed by dividing the distance between two positions 
by the time interval and assuming the errors at the two positions are 


equal, the standard deviation in wind speed due to GMD-1 errors (Oc) 


oe 


TABLE VIII 


Standard Deviation of Fosition Error (meters) 
Due to GMD-1 Measurement Errors 





can be expressed by 


eae Kee (9) 


Considering the worst case in Table VIII, the standard deviation in 
wind speed due to GMD-1 errors at 15 km would be 2.7 m sec | for a 
two minute interval and 10.8 m rca for a 30 second interval. 

Up to this point the height (Z) has been considered exact, but 
this is definitely not the case. To get more realistic values for 
wind speed errors one must incorporate the errors in range due to 
errors in height of the radiosonde. Unfortunately these errors in 
height must be multiplied by the cotangent of the elevation angle to 
be converted to range errors with the result that low elevation angles 


compound the problem. Since the errors in height, and the errors 


O97» 


in range due to GMD-l1 errors in elevation angle, Opp? can be assumed 


to be independent, the total error in range Oo, is given by 


R 


Op = op)” of (oy corey (10) 
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Using values of o, interpolated from Table VII and Opp from Table VIII, 


Z 


Op is computed according to equation (10). Values of Oy are then 


combined with Op to get Oy » 


a= Vion)” + (oy)°. ¢ll) 


the total standard deviation of position 


error according to 


The standard deviation in wind speed, Org» is then computed by 
af 
Ory a" with At = 2 minutes, and the results are tabulated in 
Table IX. 
TABLE IX 
Standard Deviation of Position Error and Wind Speed 
Due to GMD-l Measurement Errors 
and Radiosonde Height Errors (meters) 
Oy OW 
Mean Wind | Height Oo, 0, cote | Opp Op On (mervene dal) Ort ie 


ea SS Bae ae ee 
ee Le ay eg Beas [ee he TES 


Coma] = fe Pelrof | 
vanefom foe] ms [ue few psof 2 | 23 
ome] oe [ee foofaef [a 


It is obvious that radiosonde height errors when converted to 


10 m sec 


range errors are predominant in determining wind speed and position 
errors. This kind of error is minimized by making use of a precision 
radar to track balloons for upper-wind data. If a FPS-16 were used, all 
the values in Table VIII would be one-fifth as great. The five-meter 


standard deviation in slant range must be multiplied by the cosine of 
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the elevation angle to yield errors in range or displacement of about 
four meters which may be compared to values of Op in Table IX. 
Obviously, if a precision radar is available, it should be used when 
detailed structure of the wind field is desired. 

If one assumes the instantaneous wind speed at a particular level 
to be 1.5 times the mean wind, then reported wind speeds in the upper 
troposphere are often in error by + 15% and errors in the stratosphere 
may be in error by + 30% or more, This conclusion is based upon the 
fact that the two mean winds outlined above certainly do not repre- 
sent extreme conditions but rather conditions which are fairly typical 
in mid-latitudes. 

The problem of ascertaining the RMS error in direction remains. 

If one assumes that the wind affecting the radiosonde is along the 
azimuth of the instrument, then the error in direction is a result of 
only the error in azimuth angle. Given a synoptic situation with a 
mean wind speed of 10 m seem up to 300 mb (10 km) and a wind speed of 
an — at that level, then the radiosonde would be displaced 1800 
meters in two minutes; and the error normal to the wind would be 

On 2 or 27.4 meters. Thus the error in direction would be equal to 
arctan 27.4/1800 or about 0.9 degrees. However, if the wind at 300 
mb were perpendicular to the azimuth of the instrument, the error 
normal to the wind would be ame? or 118 meters causing an error in 
direction of arctan 118/1800 or about 3.3 degrees. The true wind is 
usually at some angle between these two extremes, and, for the case 
cited, the error would be about two degrees with a 30-degree angle 


between the measured wind direction and the line of sight. If winds 


were weaker or meaSured at a higher altitude, the direction error 


would increase. The converse is also true. In normal practice the 
wind direction and speed are read off a plotting board. Those errors 
have not been considered here although the total error in direction 
may be as large as five degrees with a weak wind field. Until 1968 
the wind direction was reported or coded for teletype transmission to 
the nearest 10 degrees so that it may usually be considered accurate 
to + 5 degrees. Presently, the radiosonde code provides for reporting 


the wind direction to the nearest five degrees. 
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IV. USEFULNESS OF VARIOUS ANALYSES AND COMPUTATIONS ON A MESOSCALE 


Reduction of the rawinsonde data gathered in the NSSL Oklahoma meso- 
network was accomplished by computer. The program used was essentially 
the same as briefly described by Kreitzberg and Brockman (1966). The 
pressure data were automatically smoothed by a (1, 2, 1) weighting of 
the log of pressure according to InP(m) = | (inP(m-1) + 21nP(m) + 
LnP(mt1) \/4. The heights of constant-pressure surfaces in 50 mb incre- 
ments from 950 mb were computed hydrostatically, but the heights of 
contacts, other isobaric levels, and isentropic levels were all com- 
puted by interpolation from the 50 mb incremented levels assuming that 
the temperature varied with puggee or that the lapse rate was dry 
adiabatic. Computer processing would induce less computational error 
than the graphical techniques used on the synoptic scale, and the com- 
puter output certainly provides more data than synoptic reports. 

In analyzing mesoscale features one does not expect the gradients 
to be linear. However, in the absence of known trends or data to the 
contrary it is difficult to justify a nonlinear analysis. The amount 
of smoothing that is commonly done on a synoptic scale map to account 
for "bad" data is not possible on a mesoscale because he features to 
be analyzed may be obliterated. Therefore the initial analyses are 
the result of linear interpolations of the data assuming all the data 
are correct. Some revised analyses are also discussed in the sub- 
sections to follow. 

At the time of the soundings considered, 1/700Z on 28 May 1967, there 


were no thunderstorms observed in the NSSL mesonetwork, and none developed 
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until after the soundings were terminated. A weak surface stationary 
front was shown in the northwest corner of the area at 1800Z by 
Fankhauser (1969) but there is no indication of a front at 850 mb at 
this map time. 


The station model used on the maps is: 


T H 
WYZ 
cs a 
(ddd/ff) 


where temperature (T) and dew point (T,) are reported to the nearest 

0.1C; XYZ is station identification; H is height to the nearest meter; 
em ae a) > , 

m is mixing ratio in gm kgm '; and ddd/ff is direction to the nearest 


degree and speed to the nearest meter per second. 4 


A. HEIGHT ANALYSIS 

For synoptic scale data the heights of standard isobaric surfaces 
are normally analyzed with 30-, 60-, and 120-meter differences between 
isohypes with contour interval increasing with altitude. The 850-mb 
surface is usually analyzed with 30-meter increments while pressure 
levels above 850 mb and below 100 mb are analyzed at 60-meter intervals. 
If one assumes standard deviations of height differences between radio- 
sondes to be 3, 10, and 20 meters at 850, 500, and 300 mb respectively, 
then the standard deviation of height differences range from one tenth 
to one third of the isohypse spacing and apparent inaccuracies can 
easily be smoothed. However, for the mesoscale case, the heights must 
be analyzed with a 10-meter difference between isohyses to get adequate 
definition. At 850 mb this should not be too bad, but at 300 mb the 


expected RMS error in height is twice as large as the isohypse interval. 
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In Figures 2a, b, and c the height analysis was done with the primary 
aim of drawing to the height data. Of secondary importance was making 
the contours fit the wind field. As a result of drawing to reported 
heights, the wind over at least part of all three maps may be described 
as antibaric. Furthermore, as one proceeds from the 850-mb through 
500-mb to the 300-mb maps, the antibaric flow becomes more pronounced. 
The observed wind field does not justify any closed centers on any of 
these maps; but, to provide any resolution, one is obliged to draw 
contours at about 10-meter intervals, and the height data forces the 
closed centers at 500 mb and 300 mb. 

The maximum height difference between stations on the 850-mb map 
is 21 meters; at 50 mb, 25 meters; and at 300 mb, 34 meters. I£f One 


considers the geostrophic wind equation 


and the height gradient between RIN and FSI on each of the three maps, 
then the geostrophic component of the wind normal to a line connecting 
these two stations weuldebesm29qe33% andesiiem — ake 850, 500:;eande 300 
mb, respectively. Obviously the wind is not geostrophic as there is at 
least some curvature to the streamlines; and, at 300 mb, a definite 
trough is evident near the middle of the map. Furthermore, one would 
expect some accelerations to be present in an area where thunderstorms 
develop a few hours after map time. The height differences between FSI 
and RIN are clearly excessive; and, from the analyses, one would suspect 
that FSI was reporting heights which were too low while RIN was report- 
ing heights which were slightly too high. Similar comparisons may be 


made with other stationss« 
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Figure 2a: 850 mb Height Analysis, 1700Z, 28 May 1967. 


47 


99°w : 98°W 97°W 







6°N 
~I.7 5836 5830 : 
22.9) WAT | 
5830 
2231/8 S 
589 
35°N 
-11.8 &83 
PVY 
j0.6 ':8} 
24%2/8 
5SEBHO 
7 5845 
ORIN 
0.7 
240/13 34°N 


Cheeni Scale KM 
oO $90 


Figure 2b: 500 mb Height Analysis, 1700Z, 28 May 1967. 
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Figure 2c: 300 mb Height Analysis, 1700Z, 28 May 1967. 
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Even though the wind cannot be described as geostrophic on these 
maps, the geostrophic assumption does provide a starting point for 
analysis. Accordingly, the contour spacing for adjusted height 
analyses at 850 and 300 mb was computed from equation (12) using 
observed winds for ae and a five-meter interval between contours. 
Streamlines were used as a guide in placing the contours. These 
adjusted analyses are Figures 2d and 2e, At 850 mb, differences in 
height between the two analyses are 10 meters or less while at 300 mb 
up to 20 meters difference is found. One would expect on Table VI 
that RMS differences in height at 850 mb would be about three meters; 
but from this 1700Z data, it appears that OL should be somewhere around 
five or six meters. Of course, there is a 5% chance that the differ- 
ence will exceed 20,. At higher levels, standard deviations in heights 


h 


given in Table VI appear to be accurate. At 300 mb Ca 19.1 meters, 
and this is compatible with observed apparent errors. 

Synoptic-scale height analysis of isobaric surfaces in mid-latitudes 
is a means of defining the pressure systems and the wind field. With 
the dense data network available over North America and with the RMS 
relative errors less than half the contour interval, it is relatively 
easy to obtain a good analysis and adequately portray the synoptic 
situation. However, for this mesoscale case, the amen analysis adds 
nothing to the understanding of the situation. At 500 mb the expected 
RMS relative error in reported height is the same magnitude as the 
necessary contour interval, and the situation worsens with increasing 
altitude. To assume that the winds are parallel to the height contours 
is a completely erroneous assumption. In fact, the height analysis in 


the middle and upper troposphere does nothing but confuse the situation. 


Streamlines or trajectories would be much more helpful. 


50 


9° w eS W/ IT Wi 





36°N 
© 1537 7 
WAT oy, 
4.2/ 8&7 
1540 
. 21e/u/' 
— 596 
6.4 1842 6 
COR / ie 10: 
a Wf 219/7 
/ 
6/10 
2156/1 a5 
/ 
16.3. 1532 / 15.9 15xB 
A FS} 4.0 6 
Aa q / a / 
232/49 J FR nese 
6.4 1553 | 
I5S4O Onn IS45 |. 
‘4 &.§ 6.8 | 
2 p ’ 34°N 
e2/ st? 1817 | 
19476 


Cele Scale KM 
° 50 


Figure 2d: Adjusted 850 mb Height Analysis, 17002, 
28 May 1967. Isohypse spacing is based upon observed 
winds and the geostrophic assumption. 
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Figure 2e: Adjusted 300 mb Height Analysis, 17002, 
28 May 1967. Isohypse spacing is based upon observed 
winds and the geostrophic assumption. 
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B. TEMPERATURE ANALYSIS 

The maximum temperature difference across the network for the 1/700Z 
28 May 1967 data is 1.8C at 850 mb, 3.4C at 500 mb, and 1.5C at 300 mb 
corresponding to 2.60, 4.90 and 2.1lo respectively. Here o is assumed 
to be + 0.7C om a constant-pressure surface. As.shown in Figure 3a, the 
850-mb isotherm pattern indicates a cool area averaging about 15.6C from 
southwest to northeast across the middle of the map. These isotherms 
are nearly aligned with the wind field. At 300 mb (Figure 3c) four out 
of five of these stations are enclosed by the -38C isotherm, and much of 
the thermal gradient is oriented with the wind. In contrast to the 850- 
and 300-mb levels where temperature differences were rather slight, the 
500-mb isotherms in Figure 3b reveal the large thermal gradient between 
FSI and CHK of 2.9C in 60 kilometers. How much of this difference is 
real and how much relative error is open to question. If one accepts 
o to be 0.7C, then this difference between adjacent stations would be 
4.lo. Therefore, most of the thermal gradient is either real or a 
result of reading the wrong ordinate value into the raw data. The 
latter would be a result of a gross error and not considered a part of 
relative error. There were no observed clouds extending to the 500-mb 
level at that time and none of the radiosondes encountered clouds or 
precipitation on their ascent. Therefore, radiation and conduction 
errors of the two thermistors should have been equal and not been the 
cause of this tight gradient. If one neglects the temperature report 
at CHK, a smaller and more realistic thermal gradient results (see 
Figure 3d). 

There is little qualitative use made of a temperature analysis on 


this scale. One may notice, however, that the wind does veer considerably 
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Figure 3a: 850 mb Temperature Analysis, 1700Z, 
28 May 1967. | 
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Figure 3b: 500 mb Temperature Analysis, 1700Z, 
28 May 1967. 
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Figure 3c: 300 mb Temperature Analysis, 1/00Z, 
28 May 1967: 
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Figure 3d: Adjusted 500 mb Temperature Analysis, 1700Z, 
28 May 1967. This analysis neglects the reported 
temperature at station CHK. Compare with Figure 3b. 
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between 500 mb and 300 mb, and in particular the wind change at FSI shows 
the effect of warm advection. There are two computer outputs in the 
single-sounding program which might at first glance use temperature 
data. Static energy will be considered in a latter section. The ad- 
vective temperature change is not computed by using observed temperatures 
but by use of observed winds and the thermal wind relationship assuming 
that the variation of observed wind with height is the same as the 
variation of geostrophic wind. Kreitzberg and peoetenan (1966) indicate 
that good results were obtained by this method in Project Stormy Spring. 
Most advective temperature changes (ATC) as computed from the thermal 
wind relationship and at all layers and stations were less than 0.5C per 
hour. If one compares the winds, temperatures, and distance between 
FSI and CHK at 500 mb, a much larger figure would be expected. From the 
definition of positive temperature advection, 


Bugee=)~ Weve, (13) 


the result, using a wind speed of seven meters per second, is +1.22C 
per hour. Equation (13) may be written in natural coordinates as 


ATC = - V ot. Differentiating for error analysis yields: 


d(ATC) = - d(V) ct Vv ae (14) 


where d(ATC) may be equated to RMS error in temperature advection as 
computed from the values at an isobaric surface; the first term on the 
right side, the error due to wind speed and direction; and the second, 
the error due to temperature. For example, disregard the reported 
temperature at CHK and consider only the gradient shown in Figure 3d 
where the average distance between isotherms is about 55 km. If one 


assumes the RMS error in wind speed to be one meter per second, Orp 


ao 


to be 0./C, and the average wind from FSI to TIK to be seven meters per 
second, then the expected RMS error in temperature advection due to error 
in wind speed would be 0.065C per hour while that due to error in tempera- 
ture would be 0.32C per hour. Combining the two results, assuming that 
the errors are independent of each other, would yield a total error in 
ATC of 0.33C per hour. The total expected RMS error for this example is 
mostly caused by the error in temperature and is the same order of magni- 
tude as the results using the computer program. The maximum ATC for 


this area of +0.25C per hour occurred at TIK in the 500-450 mb layer. 


C. MOISTURE ANALYSIS 

Measurements of moisture content are essential to the understanding 
and prediction of convective storms. While the analysis of moisture 
parameters at standard isobaric surfaces is of little use in itself, 
mixing ratio is often used in computing static and total energy, 
diabatic heating, and stability, to name a few. Water vapor is un- 
doubtedly the most variable of all atmospheric constituents, and its 
measurement is the least accurate. To investigate how well the water 
vapor content can be defined on this scale, analyses of moisture para- 
meters were attempted at 850, 500, and 300 mb and are listed as Figure 
4a, 4b, and 4c respectively. 

At 850 mb the mixing ratio analysis provides a pattern in which one 
may have confidence. The assumed standard deviation in reported mixing 
ratio for this level is about 0.6 gm ee Dew points for this level 
were also analyzed, but the error in dew point varies more from station 
to station and depends not only upon temperature but upon relative 


humidity as well. Although the RMS error in dew point does not vary 
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Figure 4a: 850 mb Dew-point and Mixing Ratio Analysis, 
1700Z, 28 May 1967. . 
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Figure 4b: 500 mb Dew-point Analysis, 1700Z, 28 May 1967. 
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Figure 4c: 300 mb Dew-point Analysis, 1700Z, 28 May 1967. 
Numbers under the dew-points at each station are expected 
RMS errors in dew-point. 
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mversely with the realtive humidity as does the RMS error im reported 
relative humidity, it does vary at about the same rate. That is, if 
relative humidity decreases at a constant pressure and temperature, the 
spread (in degrees centigrade) represented by + 0.1m. increases. The 
maximum difference of dew point values at 850 mb across the network was 
6.7C while that of mixing ratio was 3.8 gem eo 

By 500 mb the temperature was below OC so that expected RMS relative 
errors of + 10% relative humidity applies. This translates to about 
0.3 %m fn in mixing ratio corresponding to m= 3.0 gm a for the 
average temperature found at this level. Althoughisopleths of mixing 
ratio could be drawn, better definition is obtained by analyzing dew 
points. Therefore only dew point isopleths at two-degree intervals 
were drawn, but they must be interpreted with caution. The accuracy 
of the dew points so analyzed vary considerably and are not related to 
an arbitrarily chosen isopleth interval. Whether comparing mixing 
ratios or dew points, it is obvious that the eastern edge of the area 
contains more water vapor than the south-western corner. The expected 
RMS relative error in mixing ratio remains essentially the same over the 
entire map and is equal to 10% of m. which is dependent upon temperature. 
However, the expected standard deviation in dew point is also dependent 
upen relative humidity and varied from about + 8C at LTS to # 2C wat PVY 
with relative humidity at these stations 15% and 60% respectively. 

One glance at Figure 4c will indicate the necessity for analyzing 
dew points. The mixing ratio was reported as 0.1 gm ae with the 
exception of station WAT which reported 0.2 gm erie Relative humidity 


varied from about 15% to 35%. The resultant dew point RMS errors range 
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from + 8C at LTS to + 2.5C at WAT and are indicated at each station to 
the nearest half degree in Figure 4c. 

Relatively accurate depiction of water vapor content on this scale 
is possible at lower levels with temperatures above OC. In the middle 
and upper troposphere, errors are likely to be much larger although a 
qualitative judgement of the moisture content at the normal limiting 
temperature (-40C) seems possible. Relative humidity and dew point 
are available as computer output up to the limit of the sounding (100 
mb) for this network, but standard practice in the synoptic network 
limits information on water vapor content to dew points at standard 


levels where temperature is greater than -40C. 


D. MESOSCALE WIND ACCURACY 

Before discussing computations derived from analysis of the wind 
field some procedures used at NSSL in measurement and computation of 
the wind should be mentioned. Azimuth and elevation angles were 
recorded at 30-second intervals. There was a paca aert ty of recording 
these angles at six-second intervals, but this procedure was not used 
unless elevation angles were less than 12 degrees (Fankhauser, 1969). 
On 28 May 1967 the winds were rather light and would result in eleva- 
tion angles on the order of 25 to 40 degrees. Kreitzberg and Brockman 
(1966) indicated that the data taken at 30-second intervals were not 
smoothed but three-point interpolation was used to yield azimuth and 
elevation angles at each contact base. The detailed wind at each con- 
tact over the curved earth was computed by centered differencing of the 
values at the contacts at either side which provided about a one-minute 
average. Then the mean winds, which were printed by the computer, were 


computed by a (1, 2, 1) weighted average of the detailed winds. 
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Since the winds at any particular point are the result of averaging 
over five or six elevation and azimuth angles, some of the random errors 
in measurement of GMD-1 angles should be averaged out. The accuracy 
of the reported winds should approach that which Danielsen and Duquet 


(1966) reported (about one m soc and two degrees for light winds). 


E. ENERGY COMPUTATIONS 

Another computer output in the data supplied by NSSL is static 
energy at constant pressure levels which is defined as: 
SE = ce + gz + Lm. Total energy is also used and consists of static 
energy plus a kinetic energy term. Therefore, total energy is defined 


as: 
ae el ee Om (15) 


or 
oma 
TE = 1.004T + 0.00982 + 2500m + > x 10 (15a) 


where T is in degrees Kelvin, z is in meters, m is mixing ratio in 
gm/gm, V is wind speed in m te and TE is total energy in j em 
At 500 mb (Figure 2b) average values of temperature, height, mixing 
ratio, and wind speed substituted into equation (15) would yield: 


v)- 
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It 


TE x Tor 





1M0Ue (263) + One0oe (5880) + 2500" (OL 001) + ( 


Boe = 37 +12.) + WPO25. 
Thus one can see that the terms in the equation are arranged in 
descending order of contribution to total energy. 

The RMS error in total energy is composed of errors in all four 
parameters. The standard deviation in temperature is assumed to be 
0.7C; height (at 500 mb), 10.6 meters; mixing ratio (at 500 mb), 0.0003; 


and wind speed, 1.0m soe Therefore, the RMS error in the temperature 
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term is approximately 0.7; height term, 0.1; moisture term, 0.75; and 
kinetic energy term, 0.0005; all expressed in units of j sit The 
largest errors are contained in the temperature and moisture terms and 
are about equal in magnitude at 500 mb. However, the error in the 
moisture term at 850 mb would be twice as iarge as at 500 mb because 

: = el : 
Deere 12.0 gm kgm =_— and _ therefore Cee 0.6 gm kgm ~. The error in 
temperature is assumed to remain constant. If the errors in these four 
parameters are assumed independent of each other, the RMS error in total 
energy is given by: 


Ong = fK0.7)° + (0.1)° +°(0.75)- + (0-0005)- = 1.03 4 ae 


a 


The total energy at 500 mb ranged from 322.8 j aan at es 324.4 

j oui at COR. The maximum difference of 1.6 j apa may be compared 

to the expected RMS error of 1.03 j ei The maximum change in total 
energy across the 500-mb map is so small relative to the RMS error that 
comparisons between stations are of questionable use in other than a 
qualitative way. At 850 mb, however, the maximum difference of 8.3 


j gm is about five times as large as the expected RMS error for that 


level, and analysis at that level could be more revealing. 


F, ISENTROPIC TRAJECTORIES 

The use of isentropic trajectories for the study of vertical 
motions, moisture sources, and accelerations, to name a few, has been 
widely accepted as a valuable synoptic-scale tool. The ree kinematic 
isentropic trajectories derived from NSSL mesonetwork aera to study the 
meteorological properties of air parcels has been described by Van 
Sickle (1969). He adapted a computer program originated by MahIlman 


(1968) for use on this mesoscale to compute and draw these trajectories. 
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Nine other parameters were also computed at both ends of each trajectory 
but most will not be discussed here. 

For synoptic scale work the Montgomery stream function (M) is 
invaluable in adjusting the trajectory at each time step. However, 
Van Sickle found that in many cases isopleths of M were perpendicular 
to the wind; and, therefore, he had to compute trajectories kinematic- 
ally. Computer output of Montgomery stream function contains virtual 
temperature instead of air temperature; because the 9 values were com- 
puted using virtual temperature, and the isentropic surfaces were 
really surfaces of constant potential virtual temperature. Thus the 


equations for M and 9 become: 


0.286 
” 1000 
a=1(S z (16) 


and 


M=c T + gz = 1.004T + 0.00982 (17) 
pv V 


where Ty is in degrees Kelvin, z is height in geopotential meters above 
‘ates sal ' 
mean sea level, and Mis in j gm. Errors in temperature, pressure, 
and mixing ratio cause errors in both Q, and M. It is safe to assume 
that these errors don't completely cancel each other so that there is 
\ 
Some error in M on an isentropic surface. To make the values of M 
dynamically consistent with the winds and accelerations, Van Sickle had 
to make some rather large modifications of M relative to the total range 
of values across the mesoscale area. For example, on the 9 = 320 sur- 
face which was at 500 + 20 mb the values of M ranged from 320.70 to 
, -l -l ae 
Fl 0 Seam for a total spread of 0.35 j gm . Yet modifications of 


= ; 
up to 0.25 j gm or five times the isopleth interval were necessary to 


make the values of M consistent with observed winds and accelerations. 
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These modifications are equivalent to errors in temperature of 0.25C or 
25 meters in height which represent 0. 366,, and 0.80, respectively. 

Van Sickle (1969) also found that moisture content was a crucial 
factor in determining total energy on an isentropic surface. The contour 
patterns of total energy in his study closely resembled the mixing ratio 
field which indicates that the third term (Lm) dominates in both total 
and static energy. To change total energy by 1.0 j a requires either 
a change in temperature of 1.0C, a change in height of 100 meters, a 
change in wind speed of 45 meters m sec, ora change in mixing ratio 
of 0.4 gm Lees Obviously mixing ratio is the most sensitive and vari- 
able parameter. One could not expect either energy calculation to yield 
information which could be relied upon in other than a qualitative way. 

Kinematic trajectories on isentropic surfaces have proved to be 
valuable in studying local storms. Vertical motions obtained by Van 
Sickle were consistent with the observed growth and decay of thunder- 
storms which occurred the afternoon of 28 May 1967 in the Oklahoma 
mesonetwork. Convergence and divergence of trajectories were consistent 
with storm location. Height and temperature inaccuracies in the Montgo- 
mery stream function rendered it useless for studies on this scale, and 
computations involving moisture content did not prove to be reliable. 
Analyses and computations based upon adequate description in time and 


Space of the wind field were shown to be the most useful in this study. 


G. VORTICITY AND DIVERGENCE 
One of the most useful tools in predicting cloudiness and precipi- 
tation is the knowledge of vorticity and divergence patterns. Several 


methods of computing vertical motion are based upon these quantities. 
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Coleman (1969) attempted to adapt these synoptic-scale techniques to 
the study of a squall line in the Oklahoma mesonetwork which occurred 
on 30 May 1967. From streamline and isotach analyses at each 50 mb 
interval in the troposphere he computed vorticity, divergence, and 
vertical motions. 

Both divergence and vorticity contain terms involving spatial 
gradients of velocity. If wind speed is accurate to within one m sec. 
and the average distance between stations is 85 km, then the RMS error 


in these terms are on the order of 10 ane Coleman found diver- 


A 


gences from er to 10 soo and relative vorticities from 10. 


sec up to a maximum of 50 x one eer near the squall line. In 


some areas the computed divergences and vorticities were of the same 
order of magnitude as the RMS error,but near the squall line these com- 
puted quantities were at least one order of magnitude larger than the 
error. 

The vertical motions were obtained by two different schemes. One 
was based upon a vorticity relationship while the other was based upon 
divergence. The vertical motion computer program based upon vorticity 
did not converge because the advection term was too large due to very 
tight spatial gradients. This program had been successfully utilized 
for synoptic-scale studies. The divergence-derived vertical motion 
program revealed maximum vertical velocities of about 100 mb a 

It is interesting to note that Coleman did not intentionally smooth 
the wind data in the initial analyses. The winds on 30 May for the area 
were approximately twice the magnitude of winds on 28 May and averaged 


about 20 meters per second at 500 mb and 30 m pee near the tropopause. 


There was some averaging done after the initial computation of vorticity 
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and divergence. This smoothing amounted to giving the grid point in 
question on the pressure surface 40% weight and the four surrounding 
grid points 15% each. Vertical smoothing was accomplished by a 

(1, 2, 1) weighted average. The reported winds, of course, represent 
a weighted average over about 600 meters. 

This appears to be another method which can be successfully used 
for further understanding of severe local storms. Its validity is 
dependent upon the accuracy of the wind data and streamline and isotach 
analyses. The winds seem to provide the most reliable and useful 
information at constant-pressure surfaces for studies of this scale of 


motion. 
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V. CONCLUSIONS 


RMS errors between radiosondes assumed for both measured and 
derived parameters seem to have been borne out by several meteorologists 
who have analyzed the NSSL data. The minimum RMS relative error in 
height of standard pressure surfaces is probably about five meters but 
those derived for the middle and upper troposphere seem to be consistent 
with results. One would expect that radiosonde data from the synoptic 
network would contain slightly larger relative errors due to increased 
variability in the exposure of the instruments and the use of graphical 
techniques for data reduction. 

Wind accuracy on the synoptic scale is highly variable as outlined 
in Section IIIE. Errors in balloon position and therefore computed 
winds are directly correlated with wind speed, altitude, slant range, 
and elevation angle. By computer processing and successive averaging 
of GMD-1 angles and detailed winds the mesoscale wind ste appears to 
contain much less error than one would expect from synoptic-scale 
experience. 

Rawinsonde-derived data from the synoptic network are adequate to 
portray cyclones and other meteorological events of similar size. This 
is possible because the RMS relative errors are small relative to the 
magnitude of the meteorological changes encountered, and the data density 
in the continental United States is adequate for the task. Height fields, 
thermal advection, stability indices, and vorticity patterns are primary 
tools for the forecaster. However, on a mesoscale the RMS relative 


errors in height, temperature, and moisture parameters are too large 


oa 


relative to maximum differences observed across the mesoscale network 

to be useful. The amount of smoothing necessary to achieve reliable 
analyses appears to be the same order of magnitude as the respective 

RMS errors. Comparison of original and adjusted height fields confirms 
this although there is no direct comparison possible for other parameters. 

The wind data appears to contain the most useful information. The 
problem in isotherm analysis and thermal advection has been surmounted 
by computing advective temperature change utilizing observed winds and 
the thermal wind relationship. Streamlines and isotachs may be sub- 
stituted for the more conventional height analysis for depiction of the 
wind field on isobaric surfaces. Kinematic trajectories on isentropic 
surfaces over short (1.5 hours) periods of time are another means of 
describing the wind field. Both methods lead to computation of 
vorticity, divergence, and vertical motion all of which are essential 
to the understanding of thunderstorm dynamics. 

It would be most helpful if the measurement of water vapor content 
in the atmosphere were accurate enough to rely upon for quantitative 
analyses and computations of energies, diabatic heating, and the like. 
Some of the inherent errors in mixing ratio may be eliminated by using 
values which are averaged over 50 mb layers (cf. Frankhauser, 1969). 
One may attempt calculations of total moisture inflow into a squall 
line by advecting these mean values with the mean wind in each 50 mb 
layer. Nevertheless, there will always be considerable error in reported 
moisture parameters as long as either the lithium chloride or the carbon 
element are used. The former element is reluctant to measure high 


relative humidities while the latter will not measure low relative 
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humidities. An element which is more sensitive with an RMS error of 
<i. _. ae. ’ 
about 0.1 gm kgm in mixing ratio is badly needed for both routine 


radiosonde use and detailed mesoscale investigations. 
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VI. SUGGESTIONS FOR IMPROVEMENT IN RAWINSONDE ACCURACY 


The most important improvement to be made in radiosonde accuracy is 
to get a better humidity element in routine use. Large scale air-sea 
interaction studies would be possible and diabatic heating could be 
studied with some accuracy over land. A better sensor is even more 
important for mesoscale studies. For several years barium flouride 
and aluminum oxide elements have been under development by the National 
Bureau of Standards and Air Force Cambridge Research Laboratories, 
respectively. Both featured fast response times, and early reports were 
encouraging (cf. Jones (1963); Chleck (1966)). However, Brousaides 
(1968) has concluded that the aluminum oxide element is not suitable 
for radiosonde use because of large hysteresis, small resistance change 
with relative and specific humidity, and temperature dependence. The 
barium flouride element has been found to lack the desired stability in 
storage. It might be worth the cost and effort if a batch were manu- 
factured, calibrated, and used as an experiment in the NSSL network. 
Failing that, an accurate aircraft borne sensor could be employed for 
mesoscale research. 

The next improvement in rawinsonde data in order of priority should 
be observed winds. It is obvious that if precision radars, such as the 
 EPS-16, are available for balloon tracking, they should be used. How- 
ever, it is unlikely that such radars will be available in sufficient 
quantity in the next decade or so. The next best approach would be the 
use of a transponder-radiosonde such as that developed by the Weather 


Bureau which is capable of being used with the new GMD-2 tracking 
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equipment. Either way, position errors in range are always less than 
errors in slant range; whereas, with the present system, errors in 
height of the radiosonde are magnified in range by multiplying by the 
cotangent of the elevation angle. 

Errors in temperature and height would be next in line for improve- 
ment. A shielded bead thermistor could be suspended from an outrigger 
attached to the radiosonde as a replacement for the presently used rod 
thermistors. This would decrease solar radiation and conduction errors 
as well as reducing response time. These errors are particularly 
important in the stratosphere, but it is felt that relative errors could 
also be reduced in the troposphere with corresponding decreases in RMS 
errors in height of isobaric surfaces if data reduction were accomplished 
by computer. Smaller decreases in these height errors would be realized 
if the graphical technique were used. Camp and Caplan (1969) investi- 
gated temperature errors of a bead thermistor which was mounted on a 
JIMS PHERE and compared with standard radiosonde measurements. The bead 
thermistor showed faster response and more detail throughout each flight 
and considerably less total error at altitudes above 11 km. 

If angular measurements and slant range of a transponder-radiosonde 
system were sufficiently accurate (equivalent to the FPS-16), it would 
be practical to compute height using range and elevation angle. Height 
analysis for mesoscale studies would again be feasible. c addition, 
pressure could be computed from height and temperature data more 
accurately than it is now measured with the baroswitch. These methods 
of obtaining heights and pressures would be particularly valuable at 


altitudes above 700 mb. 
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